Mean field and pairing properties in the crust of neutron stars 
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Properties of the matter in the inner crust of a neutron star are investigated in a Hartree-Fock 
plus BCS approximation employing schematic effective forces of the type of the Skyrme forces. 
Special attention is paid to differences between a homogenous and inhomogeneous description of 
the matter distribution. For that purpose self-consistent Hartree Fock calculations are performed 
in a spherical Wigner-Seitz cell. The results are compared to predictions of corresponding Thomas 
Fermi calculations. The influence of the shell structure on the formation of pairing correlations in 
inhomogeneous matter are discussed. 



I. INTRODUCTION 



The determination of the equation of state (EoS) for nucleon matter is an important ingredient for various investiga- 
tions of astrophysical objects. A lot of attention has been paid to the EoS of baryonic matter at densities above the 
I> I saturation density of nuclear matter {po « 0.16 fm~^). These densities should be relevant to describe the interior of 
\l neutron stars and are of particular interest as they give rise for speculations about exotic phases of baryonic matter 
including the existence of kaon condensates or quark matter 

However, it is not only this regime of very high densities also the crust of neutron stars should be a very intriguing 
' phase of baryonic matter. At those lower densities and low temperatures the free energy should be reduced by allowing 
, for a phase of inhomogeneous baryonic matter. This nonuniform matter should consist out of a lattice of quasi-nuclei 
• embedded in a gas of electrons and possibly also a sea of neutrons. Thomas Fermi calculations, which are based on 
non-relativistic or relativistic mean field calculations for the homogenous system, exhibit these features 0,0, 0,0, 111- 
The presence of neutron superfluidity in the crust of neutron stars seems to be well established. This is based on 
investigations solving the BCS gap equation in homogenous nuclear matter IT(i ITU IT^ . The precise knowledge of 
O ; the pairing gap is very important for the determination of the cooling rate of neutron stars. The superfluidity of 
the material of the crust should have a significant influence on the rotation frequencies of the star. The formation of 
glitches may be related to the vortex pinning of the superfluid phase in the inhomogeneous crust. 
The aim of the investigations presented in this paper is to explore basic properties of this inhomogeneous phase of 
baryonic matter on a level beyond the Thomas Fermi approximation. What are the consequences of shell effects in 
the quasi-nuclci in the systems for the bulk-properties of this material? How do they influence the proton fraction in 
' the P equilibrium? What are the consequences of these quasi-nuclei on the pairing properties of the system? 

To answer this kind of questions we perform nuclear structure calculations in a Wigner-Seitz cell of spherical shape. 
The assumption of such a spherical cell is not designed to allow for inhomogeneities of the matter which is different 
from the formation of quasi-nuclei or bubbles in the homogenous matter. Therefore we do not consider a formation 
of rods or slabs of increased density 0|- Also at this stage we do not aim at a complete survey of this matter covering 
a large region of various densities. Instead we wish to explore some features of nuclear structure calculations beyond 
the Thomas Fermi approach for a few selected examples. 

For that purpose we perform Hartree-Fock and Hartree-Fock plus BCS calculations in an appropriate basis of the 
spherical Wigner-Seitz cell. For the NN interaction we use rather simple effective forces of the Skyrme typelT^, Tlj. 
Special attention is paid to the pairing gaps resulting from calculations of the homogenous matter and the non- 
homogeneous system. 

After this introduction we discuss in section 2 the techniques and results of Hartree-Fock calculations for the /3-stable 
matter within a spherical Wigner-Seitz cell. The determination of pairing properties is presented in section 3, which 
is followed by a short summary and conclusions in section 4. 
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II. HARTREE-FOCK CALCULATIONS IN A SPHERICAL BOX 



The single-particle wave functions for the nuclcon considered in our calculations are expanded in a complete basis of 
orthonormal states defined within a spherical box of radius R. Such a orthonormal set of basis functions, which are 
regular at the origin in the center of the box is given by 

$»ijm(r) = {r\iljm) = Riiir)yij,ni'&, (p) ■ (1) 

In this equation yijm represent the spherical harmonics including the spin degrees of freedom by coupling the orbital 
angular momentum I with the spin to a single-particle angular momentum j. The radial wave functions Ru are given 
by the spherical Bessel functions, Ra{r) ~ ji{kir) for the discrete momenta ki, which fulfill 

RaiR) ^ NajiihR) = . (2) 

The normalization constant 

forZ > 0, 
for Z = . 



ensures that the basis functions are orthogonal and normalized within the box 







Adopting this basis of eigenstates for the kinetic energy one can try to describe homogenous nuclear matter with a 
constant density p^r and p,y for protons and neutrons, respectively, by filling all basis states with momenta ki below 
the corresponding Fermi momentum kp,^ and kp^. The local density is then given by 

pir)=Y,e{kF-k,)i2j + l)Rlir), (5) 

ilj 

where 9 stands for the Heaviside function. Results for this local density considering a spherical box with radius R—25 
fm and trying to describe a system with a density of 0.06 fm""^ is displayed in Fig.^ From this figure one can see that 
the local density is in reasonable agreement with the mean value except at the border of the cavity at r = i?. Since 
all wave functions by construction are bound to disappear at this borderline, the same is also true for the density. 
In order to cure this deficiency one could try an alternative basis for the radial functions by the boundary condition 
that the radial derivative vanishes at the surface of the box 

Ra{r) ^ NaMk^r) with — ^(i?)=0. (6) 

or 

So all these basis function exhibit an extremum at the surface of the Wigner-Seitz cell which leads to a maximum of 
the local density at ?■ = i? if these functions are employed in (O as on can see from the curve labeled "Basis 2" in 

Fig.n] 

Bonche and Vautherin|l5| suggested to use a mixed basis by employing e.g basis states with the boundary condition 
(PI for states with even / and those with the boundary condition © for states with odd orbital angular momentum /. 
This recipe cures the deviation of the local density from the mean value at the surface (see also Fig. ^ . Nevertheless 
the representation of the homogenous phase of nuclear matter within a spherical box of finite size leads to slight 
fluctuations in the local density. Furthermore one must be aware of shell effects, which are due to the finite size of 
the spherical box considered. 

In order to explore the influence of these shell effects on the evaluation of expectation values, we have considered 
spherical boxes of various radii and calculated the binding energy per nucleon for homogenous nuclear matter using 
the representation of the plane wave states discussed above. We have considered values for the radius R ranging from 
15 fm to 25 fm. Most of the results turned out to be rather insensitive to the choice of the basis (((2J), ® or mixture). 
Therefore, if not stated differently, we will present results for the basis ||2J) only. 

As a first example we consider the homogenous system of neutrons, protons and electrons in /3 equilibrium. This 
means that for any value of the baryon density considered 



(7) 
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we determine the proton abundance 



(8) 



P 



by the requirement that the Fermi energy for the neutrons is identical to the Fermi energy of the protons plus the 
Fermi energy of the electrons with the density of electrons being identical to the density of protons. 
Results for the energy per nucleon and the proton abundance are displayed in figures |21 and |31 respectively. For 
the range of densities considered the results obtained in the spherical boxes of different size agree rather well with 
the corresponding observables calculated in the infinite system. A significant discrepancy is only observed in the 
calculation of the total energy at higher densities considering boxes with small radii. In these cases the surface effects 
displayed in Fig^ as well as the shell effects lead to energies which are to small as compared to the result for the 
infinite system. 

After this test of the box basis to describe the homogenous system, we now turn to the Hartree-Fock (HF) description 
within such a spherical box. The HF single-particle wave functions are expanded in the basis of (Q) by 



The number of basis states A'' is chosen to guarantee that the results are not affected by this limitation. Results of 
self-consistent HF calculations are displayed in Fig. 0] showing the density of protons and neutrons as a function of 
the radial distance from the center of a box with radius R=21 fm. the particle numbers for protons {Z=20 in this 
example) and neutrons (A^=388) have been determined to fulfill the condition of f3 stability, employing the Fermi 
energies evaluated in the HF calculation. 

The density distribution for the protons is different from zero only near the center of the Wigner Seitz cell. In this 
region also the density of the neutrons is considerably larger than in the rest of the spherical box. Therefore we 
can interpret this configuration as a quasi-nucleus embedded in a neutron see. The single-particle energy spectrum 
for the protons exhibit a clear shell structure up to energies around the Fermi energy. The density of states around 
the corresponding Fermi energy is larger for the neutrons. In fact it is close to the one of the homogeneous matter 
described in the spherical box of this size. 

For the example of Z—2Q protons and A^=388 neutrons in a Wigner Seitz cell with radius R=21 fm considered in 
Fig. ^ we obtain a global density of p = 0.0105 fm~'^ and a abundance for the protons oiYp = 0.049. The energy per 
nucleon for this self-consistent solution of the HF equations is about 2 MeV per nucleon below the energy which was 
obtained for the homogenous distribution of matter in /3 equilibrium. 

Fig. m exhibits this gain in energy due to the formation of quasi- nuclei in /3-stable nuclear matter. For small densities 
around p = 0.01 fm^'^ the HF solution with localized quasi-nuclei yields an energy which is about 3 MeV per nucleon 
below the energy of the corresponding homogenous matter. This gain in energy decreases with increasing density to 
around 1 MeV per nucleon a,t p — 0.04 fm^"^. The results are rather independent on the size of the Wigner Seitz cell 
under consideration. 

The Hartree-Fock calculations allowing for an inhomogeneous distribution of matter, however, also yield different 
results for the proton abundances as compared to the results obtained for homogenous nuclear matter in (3 equilibrium 
(see Fig. Thomas Fermi calculation, which allow for an inhomogeneous distribution of matter, show proton 
abundances which are slightly above those derived from the calculation of homogenous matter. This indicates that 
the fluctuations in the density, which are taken into account in the Thomas Fermi calculations, tend to enhance the 
proton abundances, as regions with high densities contain a larger fraction of protons. However, the main effect in the 
enhancement of the proton abundances observed in the Hartree Fock approach is not described by the Thomas Fermi 
model. It originates from the pronounced shell structure of the proton single-particle energies. The Hartree-Fock 
calculation favors distributions of matter with a quasi-nucleus showing a closed shell for the protons. 



After we have noticed the effects of shell structure on the decomposition of the nuclear material, we would like to 
address in this section the question how this shell structure can affect the pairing properties of the nuclear material. 
In particular we consider neutron- neutron pairing for neutron pairs with total momentum equal to zero in ^Sq partial 
wave for the relative motion. Using the standard BCS approach for homogenous matter the pairing gap A(fc) for a 
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pair of neutrons with relative momentum k = |fc| is obtained by solving the gap eauation[9l IllL 

A{k)=-- r dk'k'^v{k,k') — ^^^'^ — (10) 



2^(4-£^r + A(fc')2 

Here V{k, k') denotes the matrix elements of the NN interaction in the ^Sq partial wave, £k the single-particle energy 
for a nucleon with momentum k and ep the Fermi energy. Instead of using the matrix elements for a realistic NN 
interaction which is fitted to the NN scattering data, we have chosen to employ an effective interaction which is of 
zero range, density dependent and has been proposed by Bertsch and Esbensen [TgLIitI 



V{ri,r2) = Wo 



1—7/ 



\ Pa 



5{ri~r2), (11) 



where vq, tj and a are adjustable parameters and po is the saturation density of nuclear matter. For such a zero- 
range interaction a cut-off parameter ec must be introduced in the gap equation to truncate the integral to states 
with energy Sk less than ec- These four parameters can be constrained to reproduce the TV TV scattering length and 
Garrido et al. [ig determined various sets of parameters which reproduce pairing gaps for nuclear and neutron matter 
calculated with realistic NN interactions. We used the parameters vq = 481 MeV fm^, 77 = 0.45, a = 0.47 and ec — 
60 MeV and verified that this set of parameters leeds to pairing gap s for homogenous neutron matter, which arc in 
fair agreement with those evaluated for the CD Bonn interaction 19]. 

We now turn to the solution of the HF-I-BCS equations ,2fl|, which in the spherical Wigner-Seitz cell have the form 

{Snlj — £F)Unlj +\iljVnlj = EnljU^ij 
— {Snlj — ip) Vnlj +^nljUnlj — EnljVnlj , (12) 

where 



is the energy of the quasiparticle state with quantum numbers n,l,j and e^ij the corresponding single-particle energy 
determined in the HF equations. From the coefficients Unij, v„ij and the corresponding single-particle wave functions 
"^nijm (see eq.®) one can calculate the anomalous density 

x(r)=5:(2j + l)^^^*L,.(r), (13) 

nlj 

and the expression for the normal density is changed into 

p{r)^J2'^2j + l)vli^^li^{r). (14) 

nlj 

Using an interaction of zero range like 11111 leads to a definition of a local gap-function 

A(r) = Vir)x{r) , (15) 
from which one can calculate the state-dependent pairing gaps 

A„y = I A{r)^li^{r)rUr. (16) 

The BCS equations (|12|l - H16|l have to be solved together with the HF equations in a self-consistent way. The Fermi 
energy is adjusted to reproduce the requested particle numbers or densities for protons and neutrons. These 
equations can be reduced to an approach, which we will call BCS with plane waves in the spherical box, by restricting 
the radial wave functions for the single-particle states to the spherical Bessel functions defined in 101 or ((Bj). 
Results for the gap function A(r) defined in H15|l for /3-stable matter of a density p = 0.02 fm^^ are displayed in Fig.d 
The gap function evaluated for the plane wave basis are presented by the dashed line. We find that this gap-function 
fiuctuates around 1.8 MeV with a sharp drop at the boundary of the spherical box. This drop is of course related to 
the boundary condition (0) of all wave functions. As already discussed for the local density above, this deficiency can 
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be cured by using a mixed basis, switching between ^ and lO for states with even and odd parity (see lower panel 
of Fig.ig. 

Solving the BCS plus Hartree-Fock equations one obtains a gap-function with values which are suppressed by about 
25 percent in the area of the quasi-nuclci. From this result one may speculate that the inhomogeneous matter leads to 
regions of enhanced densities in which the formation of pairing correlations is suppressed to some extent. Therefore 
these regions might be considered as nuclei for the formation of normal matter within the superfluid phase of neutron 
matter. Also one could imagine that these quasi-nuclei could lead to vortex pinning in the rotation of the superfluid 
crust of neutron stars. 

Before drawing any conclusions of this kind, we would like to compare the results for the state-dependent gaps A„ij 
displayed in Figs. |H1 and El for the plane wave plus BCS and HF plus BCS calculations, respectively. Comparing these 
results, one finds that plane wave plus BCS calculation leads to a gap Anij which is almost independent of the state. 
The calculation of the inhomogeneous system on the other hand leads to a few single-particle states, which are much 
more bound than the lowest states in the homogeneous calculation (-42 MeV as compared to -9.5 MeV). These deeply 
bound states are localized near the origin. This large binding of the single-particle states which leads to a reduction 
of the density of states is also responsible for the reduction of the pairing gap for these low-lying states. On the other 
hand, the states with single-particle energies close to the Fermi energy ep show a value for the gap parameter which 
is very close to the one derived for the homogenous system. 

This implies that the reduction of the gap- function A(r) in the region of the quasi-nuclei discussed above manifests 
itself mainly in a reduction of the pairing gap for the deeply bound single-particle states. The evaluation of response 
functions and other observables, however, is in general sensitive to the gap at the Fermi energy. Therefore one cannot 
expect major differences in the behavior of the inhomogeneous as compared to the homogenous system of matter with 
respect to the pairing properties of the neutrons. This shall be different for smaller densities, where the Fermi energy 
for the neutrons drops to values at which the shell structure of the single-particle energies is significant. 

IV. CONCLUSIONS 

The transition from homogenous to inhomogeneous matter as it occurs in the crust of neutron stars has been in- 
vestigated. Special attention has been paid to the consequences of shell effects, which occurring assuming a phase 
of quasi-nuclei embedded in a sea of neutrons. For that purpose Hartree-Fock plus BCS calculations are performed 
assuming a basis of single-particle states for a spherical Wigner-Seitz cell. It is observed that the shell effects lead to 
a substantial increase for the proton abundance in /3-stable baryonic matter. 

Shell effects are also responsible for a decrease of the localized pairing gap in the region of the quasi-nuclei. A more 
detailed analysis, however, shows that this reduction of the local gap function A(r) is mainly due to a reduction of 
the pairing gap for the deeply bound single-particle states. The pairing properties for the single-particle states close 
to the Fermi energy are similar to those obtained for a homogenous description of the system. 

The present studies are based on simple parameterizations for the effective NN interaction. It may be of interest to 
study whether the predictions for the transition from homogenous to inhomogeneous matter also hold, when more 
realistic NN interactions are employed. 

The use of a spherical Wigner-Seitz cell is a source of various deficiencies. The boundary conditions at the border of 
the cell lead to fluctuations in the density, which complicates the comparison with the infinite homogenous system. 
This could be improved by employing a Cartesian basis, which is more involved from the numerical point of view. 
The use of a Cartesian basis, however, would also allow the study of non-spherical structures. 
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FIG. 2: Energy per nucleon for homogenous matter in /3 equilibrium. Results for infinite matter are compared to those obtained 
in spherical boxes of different radii. The Skyrme I has been used for the NN interaction. 




FIG. 3: Proton abundance Yp for homogenous matter in l3 equilibrium. Further comments see Fig. |5| 
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FIG. 4: Density distributions derived from Hartree Fock calculations for nuclear matter in /3 equilibrium. 




FIG. 5: Energy per nucleon of /3-stable nuclear matter in the homogenous phase (symbols connected by dashed line) and 
inhomogeneous phase (symbols connected by solid line). 




FIG. 6: Proton abundances of /3-stable nuclear matter in the homogenous case of homogenous infinite matter (solid line), 
inhomogeneous matter determined in the Thomas Fermi approximation and in Hartree-Fock calculations employing Wigner 
Seitz cells of different radii. 
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FIG. 7: Local gap function A(r) as defined in HSU for /3-stabIe matter of density p = 0.02 fm~^. Results are displayed for the 
plane wave plus BCS approach (dashed line) and the HF plus BCS approach (solid line) While the upper panel exhibits results 
obtained in Basis 1 (see eq.Q, the lower part of the figure displays results obtained for the mixed basis. 
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FIG. 8: Results for the state-dependent pairing gaps A„;j defined in for /3-stable matter of density p — 0.02 fm~^ as 
derived from plane wave plus BCS calculation. These gaps are displayed as a function of the corresponding Snij- Each gap 
is presented by a symbol, gaps belonging to states with identical orbital angular momentum / and angular momentum j are 
connected by a line. 
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FIG. 9: Results for the state-dependent pairing gaps A„ij as derived from a HF plus BCS calculation. Further details see 
Fig.lHl 



